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ABSTRACT 



We report the discovery based on GLIMPSE data of a proto-stellar system driving a bipolar outflow . The bipolar outflow closely 
resembles the shape of an hourglass in the infrared. The total luminosity of L, OM ; = 5507 L Q , derived from IRAS fluxes, indicates 
the ongoing formation of a massive star in this region. The spectral energy distribution (SED) of the driving source is fitted with an 
online SED fitting tool, which results in a spectral index of about 1.2. This, along with the presence of a bipolar outflow, suggests the 
detection of a Class I protostar. The driving source indicates prominent infrared excesses in color-color diagrams based on archived 
2MASS and GLIMPSE data, which is in line with an early evolutionary stage of the system. 
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The mass of a young stellar object (YSO) is the key fac- 
tor that determines its evolutionary path after its formation. 
Young stellar objects with masses below and above 8 M Q 
experience substantially different proc esses in their formation 
(IZinnecker & Yorkdl2007l: lEvansll201 ll) . For low-mass star for- 
mation^ widely accepted scenario has been proposed by 
IShu et all d!987l) . in which four distinct evolutionary stages were 
defined. On the basis of the st ar formation scenario proposed by 
IShu et alj dl987h . lLadal d!987l) developed a widely used classifi- 
cation scheme for YSOs primarily based on their spectral energy 
distributions (SEDs) and the spectral index at the infrared. With 
an evolutionary sequence fro m early to late type, YSOs were 
classified into Class I to III ( lLadalfl987l) . Later on, an earlier 
evolutionary phase called Class was added to this scheme by 
lAndrg(fT99l . 

H owever, the scenario of ma ssive star formation remains un- 
clear (IZinnecker & Yorkell2007l) . Two competing schemes are 
well known in the literature. As their low-mass counterparts, 
massive stars were proposed to form via accr etion of ambient 
gas, t hough the process is much more violent dGarav & Lizand 
On the other hand, massive stars were alternatively pro- 
posed to form by m erging of low-mass stars at the center of 
dense stellar clusters dBonnell et al.lll998h . As observational ev- 
idence accumulates, the accretion scenario turns out to be the 
more likely one. 

Outflows are commo nly associat e d with both low - mass and 
massive star formation dLadal 119851: iBachillerl fl996t IWu et alj 
l2004t lArce et alj|2007l) . Massive outflows, as the feedback of 
star formation, may create strong impact on the ambient condi- 
tions such as chemical abundance and density, the final mass of 
the central star in formation and star-forming activities in both 
of the close vicinity and far reaches. 



In this paper, we report on a newly discovered Class I-type 
protostar driving an hourglass-shaped outflow. Little is known 
in the literature about the likely exciting source GLIMPSE 
G012.4013-00.4687, located at (R.A.=18 A 14 m 24 I .56, 
Dec.=-18°24'42".69 J2000). The associated 2MASS point 
source given in the highly reliable GLIMPSE catalog is 
2MASS J18142457-1824421. And the 2MASS photometric 
data reveal high infrared exce ss. An IRAS source (IRAS 



1811 4-18251 dCassaretalJll989h is located 9.3" to i ts north- 
west. ICodella et alJ dl995l) and MacLeod et alJ dl998l) tried to 
search for water and 6.7 GHz methanol masers toward IRAS 
18114-1825 but achieved negative results. The nature of IRAS 
181 14-1825 thus remains unclear. 

This paper is arranged as follows: we present a description of 
the archival data and observations in Section|2] while the detec- 
tion of the bipolar outflow system and the subsequent analysis 
are described in Section[3] In Section|4]we exclude the possibil- 
ity of an evolved object. The evolutionary stage of this source 
is investigated in Section [5] and the results are summarized in 
Section|6] 



2. Date acquisition and analysis 

2.1. Archival data 

GLIMPSE (the Galactic Legacy Infrared Mid-Plane Survey 
Extraordinaire) is a le gacy project of the Spitzer Space Telescope 
dWerner et al.l |2004|) . The Infrared Array Camera (IRAC) 
(Fazi o et a IT l2004h was employed in this survey, which covers 
four infrared wavebands centered at 3.6, 4.5, 5.8, and 8 pm, 
respectively. Archived GLIMPSE I dBeniamin et al.ll2003l) data 
were used in this work. A composite color image was compiled 
from the imaging data at 3.6 pm (blue), 4.5 pm (green) and 8 
pm (red). Photometric data from the highly reliable GLIMPSE I 
catalog were used to fit the SED of the exciting source. We only 
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considered sources with data available in all of the four IRAC 
bands. 

Archived data from the 2MASS Point Source Catalog (PSC) 
were also used in our work. To guarantee the relia bility of the 
data, w e adopted the source selection criteria from lLi & Smith! 

These include that a) each source extracted from the 
2MASS PSC must have a certain detection in all J,H & Ks 
wavebands; b) only sources with a Ks-band signal-to-noise ratio 
above 15 are selected. The retrieved 2MASS PSC data were used 
in the SED fitting and to obtain the near-infrared color-color di- 
agram. 

T he Bolocam Galactic Pl ane Survey (BGPS, iBallv et al.l 
l2010tlRosolowskv et al.ll20Toh is a 1.1 mm continuum survey of 
the Galactic Plane performed with the Bolocam on the Caltech 
Submillimeter Observatory, which is operated by Caltech under 
a contract from the NSF. A cutout image from archived Bolocam 
1 . 1 mm data was employed in this investigation. 

Photometric data from MSX, DENIS and IRAS were also 
included in the construction and the fitting of the SED of the 
central exciting source. 

2.2. Optical spectroscopy 

Low-resolution spectroscopy of the central source of the bipo- 
lar outflow was performed on the nights of 11 and 13 
September 2010 with the 2.16 m telescope of the National 
Astronomical Observatories, Chinese Academy of Sciences. An 
OMR (Optomechanics Research Inc.) spectrograph and a PI 
1340x400 CCD were used. Here, a 200 AmirT 1 grating and a 2" 
slit, across the source along the declination direction (PA.=0°), 
resulted in a two-pixel resolution of 9.6 A. A standard star, HR 
7596, was observed as the flux calibrator. Because the source 
was faint in the optical, we took three exposures, one on 11 
September and the other two on 13 September, and then com- 
bined them to get a better signal-to-noise ratio. The exposure 
time lasted an hour for each exposure. 



3. Results 

Fig-E a ) presents the Spitzer/IRAC 8 fim mosaic of a 25'x25' re- 
gion, in which the outflow structure is marked with a blue box. It 
is located at the junction of two large filaments. These filaments 
are recently reported as potential regions of active star formation 
in theoretical and observational studie s ( Dale & Bonnellll2~01 U 
lAndre et al.ll20Tot iMolinari et all2010t iHill et alJ^oflO n Fig. 
Q2 a X the eastern extended cavity with a bright central star and 
the southern filamen tary structure are parts o f a complex bubble, 
N6 as presented by iDeharveng et all (l2010l) . The nature of N6 
is still unrevealed. The bipolar source of interest here is situated 
between the two components of N6 and looks to be linked to 
both with dust showing week PAH emissions. This could be due 
to the projection effect. Additional observations with molecular 
lines at longer wavelength are necessary to make this clear. 

As shown in the composite image (Fig. []]£>), compiled with 
IRAC 8 jjm (red), 4.5 /mi (green), and 3.6 //m (blue) mosaics, the 
GLIMPSE data reveal an hourglass-shaped structure with a NE- 
SW orientation (PA. =52°). The prominent extended emission 
of 4.5 fim, containing both H2(v=0-0, s(9, 10,11)) lines and the 
CO(v=l-0) band head, suggests the detection of a massive young 
stella r object with ongoing outflow activity dCvganowski et al.l 
120081) . The probable exciting source, GLIMPSE G012.4013- 
00.4687, is marked with a magenta cross, which is located just at 
the edge of the error ellipse of IRAS 181 14-1825. This suggests 



an association of the exciting source of the bipolar structure with 
IRAS 18114-1825. 

The bipolar structure is also detected at SpitzerfMlPS 24 /urn. 
The Bolocam 1 . 1 mm dust continuum reveals an elongated struc- 
ture perpendicular to the bipolar system, which is interpreted 
here as a massive envelope and will be discussed in detail in 

3.1. Distance 

The interstellar extinction to the infrared point source is calcu- 
lated to be Ay = 2.384 mag, according to its 2MASS colors (e.g., 
J - H and H - K s ). A distance of 2. 15 Kpc is derived based on 
an extinction model tow ard the Galactic coordina te of the bipo- 
lar system developed by lAmores &Lemnel(l2005l) . 

With the extinction law (AkJAv = 0.112) proposed by 
iRieke & Lebofskvl (TT981 . we obtained the extinction in the K s 
band, Ak s - 0.267 m ag. A distance of 2.67 Kpc was estimated 
based on the model of [Marshall et al.l(l2006l) . The averaged value 
(d=2.41+0.26 Kpc) of the distance estimated with the two dif- 
ferent models was then adopted. 



3.2. Physical parameters 

The total flux density of the infrared point source c an be derived 
from the IRAS data with the following equation dCasoli et al.l 
119861) : 



F(l(n 13 Wrrr 2 ) 



1.75 



^12 . ^25 . ^60 , ^100 

0.79 2 3.9 9.9 



(1) 



where F12, F25, ^60, and F100 are the flux densities in Jy at 12 
fim, 25 fim, 60 fim, and 100 fim respectively. With the estimated 
distance of 2.41 Kpc, we obtain a total luminosity of L tota i = 
5507 + 1200 L . A luminosity this high suggests the ongoing 
formation of a massive star or a quite massive stellar system in 
this region. 

A dust temperature of T</ a 29.8 K was es t imate d from the 
color temperature T e defined bv lHenning et al.1 (Il990l) . 



T d * r f (60/100) 



96 



(3+/ j)ln(M>)_i n( ^) 



(2) 



where and F100 are the flux densities in Jy at 60 fim and 
100 fim respectively, and B the emissivity index of dust particles, 
which is set to be 2. 

A BGPS source is detected with Bolocam at 1 .1mm (see Fig. 
Q]c). With an assumption of an isothermal condition in this re- 
gion, the m ass of the BGPS source is obtained by adopting equa- 
tion (10) of lRosolowskv et alj d2010l) . 



M : 



d 2 S 

By(T)K v 



13.1M ( 



1 Kpc 



) 2 (ff)[ 



exp( 



13.0 K 



)-l 



Uy' L exp(JM)_i 



-], 0) 



where S v is the flux density measurement in the BGPS catalog, 
B v is the Planck function an d the opacity k v is assumed as 0.01 14 
cmV dEnoch et al.ll2006l) . The estimated mass of the BGPS 
source is about 69 + 15 M , which is in line with the estimated 
mass of the envelope (M em , « 22 M , see j35.21 i by the SED 
fitting. This substantiates the detection of a massive envelope in 
association with the outflow system. 
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Fig. 1. (a):Spitzer/lRAC 8 fim mosaic of the bubble N6 in lDeharveng et all d2010l) . The blue dashed box indicates the the outflow 
structure, (b): The composite color image of the IRAC mosaics of the bipolar system, where the images at 8, 4.5, and 3.6 //m 
are mapped in red, green and blue respectively, (c): Gray-scaled color-inverted Spitzer/MlPS 24 //m image is overlapped with the 
contours of Bolocam 1.1 mm continuum emission. The intensity of the contours varies from 0.1 Jy/Beam to 0.8 Jy/Beam with 10 
levels. The error ellipse of the IRAS point source (IRAS 181 14-1825) is over-plotted on both panels. The magenta cross marks the 
central source of the system. 
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4. Young stellar object vs. evolved stellar object 

We cannot exclude, however, its possible origin in a young 
planetary nebula (PN) or a proto-planetary nebula (PPN) be- 
cause these also show similar bipolar morphologies in the in- 
frared dBalick et al.|[l987HKwol[l99l . Recombination and for- 
bidden line emission from hydrogen and metals, especially H a 
and fOnil, are significant spectr al features of a PN in the optical 
dOsterbrocU 1 964t iMillerll 1 9741) . However, very weak H ff and no 
[Om] emission were detected in the optical spectrum (see Fig. 
13, which indicates that the bipolar system is not a PN. On the 
other hand, a mass estimate of about 9.2 M of the central excit- 
ing source (see 95.2b is much high er than that e xpected for a PN 
or a PPN, which is usually ~1 M (lKalerlll985l) . 

The major diameter of the bipolar structure is about 110", 
corresponding to a physical scale of 1.29 pc at the distance of 
2.41 kpc. Given an expansion velocity of 30 km s _1 (the typi- 
cal speed of exp ansion of a PPN ranges from a few to 30 km 
s^'. lKwo^l 1993b . a dynamic age of 2.096 x 10 4 yr is calculated, 
which is inconsistent with the age of a PPN that rang es from a 
few hundred to several thousand years (lKwoklll993l) . Another 
supporting evidence is the mass of the envelope derived from 
the Bolocam 1.1 mm data (see §3.2l i. which is far beyond that 
possible for material ejected from an evolved star. 

Compared with that of their pre-main-sequence counter- 
parts (Class II and III), the optical spectra l properties of pro - 
tostars (Class I) have seldom been studied d White et al.ll2007l) . 
IWhite & Hillenbrand! (12004]) carried out spectroscopic investi- 
gation into 15 Class I YSOs in the optical domain, in which 
strong emission lines (e.g., H a ) are detected. And P-cygni pro- 
files of broad H a emission of some sources reveals the features of 
outflows or jets initiated from t he cen tral protostars or disks. In 
some of I White & Hillenbrand! (l2004l) 's sample, forbidden lines 
were also detected, e.g., [S n] /16731, whose equivalent width is 
always used for the estimation of mass outflow rates. In our ob- 
servations of the likely exciting source, only weak H„ is traced, 
and no forbidden emissions (e.g.,[S n] /16731) are detected (see 
Fig. [2). The deep embedded and young nature of the massive 
protostar could be responsible for the weakness of H„ and the 
absence of [S n] A673 1 . Observations with a higher sensitivity 
and a slit along the outflow axis (PA. =52°) are needed for a 
closer investigation into its optical spectrum. 



Fig. 3. Fitted SED of the protostar. The solid line shows the best- 
fitted SED model. Variant symbols represent near- to far-infrared 
photometric data of the protostar, while the dashed line is for 
the SE D of the forthcoming s tellar photosphere in the best-fitted 
model dRobitaille et al.ll2007l) . 



5. Evolutionary status 

5.1. Implications of the morphology 

Previous studies indicate the opening angle of the outflows as- 
soci ated with massive YSOs tends to increase wit h evolutionary 
age dArce & Sargentl2006l:ISeale & Loonevl2008l) . For example, 
outflows associated with Class sources used to be highly col- 
limated, while those in associat ion with Class I protostars are 
found to have poor collimation dArce et al. I l2007h . The opening 
angle of the outflow reported here is estimated to be about 101°, 
and the collimation factor (the ratio of major to minor radii) is 
about 2. This poor collimation, together with the large opening 
angle, indicates that the central YSO cannot be in its very early 
stages of evolution. 

The investigati on of a jet-like bipolar outflow was reported 
by IVarricattl (1201 ll) . Near-infrared emissions, at JHK and IRAC 
4.5/im wavebands, reveal its well-collimated feature, which was 
interp reted as evidence for its accretion-driven nature (IVarricattl 
1201 jj). In contrast, the outflow discussed here is not well- 
collimated, and its engine is still unclear. One possible ori- 
gin of poorly collimated outflows like the one presented here 
is the entrainment in the adjacent medium from hi gh-velocity 
jets o riginated from disk winds or x-winds (e.g., lArce et all 
120071) . However, recent numerical studies suggest different driv- 
ing mechanisms for o u tflows with and witho ut good collimation 
(IMachida et al.ll2008l) . IMachida et all (l2008h found two distinct 
flows: a low-velocity flow with a wide opening angle, driven 
from the adiabatic core; and a high-velocity flow with good col- 
limation, driven from the protostar. And this scenario is sup- 
ported by other previous t heoretical works (e.g., lTomisakal2002T: 
iBaneriee & Pudritzl2006l) . And the hourg lass-shaped outflow re - 
ported here resembles that presented in IMachida et al.l d2008l) . 
Observations with outflow-tracing molecules (e.g., CO, SiO) 
will give clues to its genuine driving mechanism. 

5.2. Spectral energy distribution 

Based on the multi-wavelength photometric data, the SED of the 
central YSO was calculated and fitted with an online tool 1 de- 
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Fig. 4. Near- to far-infrared colo r-color diagrams, (a): (J - H) vs. (H - K s ) diagram. Green solid lines are th e loci of the main- 
seque nce dwarfs and giant stars dBessell & Bretj[T988h . The arrow shows a reddening vector of Ay = 5 mag dRieke & Lebofskvl 
Il985h . The dot-dashed line indicates the locus of dereddened T Tauri stars (Mever et ali i 19971) . The dashed lines, which are drawn 
parallel to the reddening vector, define the reddening band for normal field stars and T Tauri stars. Crosses are o verplotted with 
an inte rval corresponding to Ay = 5 mag. (b): [3.6]-[4.5] vs. [5.8]-[8] diagram. The classification s cheme is based Megeath et alJ 
(12004 and lAllen et all d2004 l2007h . The arrow indicates an extinction of Aks = 5 mag according to llndebetouw et al] (120051) . (c): 
log[F60/F25] vs. log[F25 /F12] diagram. Th e sample of Hn region s , T Tauri stars, Herbi g Ae/Be stars, post A GB s tars, and planetary 
nebula e are adopted from lPalla et~al1 dl99ll) . IZhang et al.l (l2006al) . IZhang et al.l (l2006bl) . ISuarez et al.l (120061) . and lTaiitsu & Tamural 
dl998l) . respectively. The protostar is presented as an asterisk in all three panels. 



veloped bv lRobitaiUe et al.l(l2006il2007h . lRobitaille et all d2006) 
presented a grid of radiation transfer models of YSOs, covering 
a wide range of possible evolutionary stages and stellar masses. 
The grid consi sts of 20,000 YSO mo dels resulting in a total of 
200,000 SEDs dRobitaille et al]|2006h . 

The best-fitted SED of the central protostar is shown as a 
solid line in Fig. [3] T he res ulting SED resembles that of a Class 
I object well dAndrelll994l) . The fitting also indicates a central 
stellar mass of ~ 9.2 ± 0.8 M Q and a total luminosity of (4 + 1) x 
10 3 L , consistent with the result of 5507 L Q estimated from the 
IRAS fluxes. The spectral index at infrared is cal culated to b e 

1.2, which is consistent with that of Class I objects lLadal fl987l) . 

5.3. Color-color diagrams 

Fig. |4] (a) is the (J - H) vs. (H - Ks) color-color diagram of 
the 2MASS point sources in a 30' x 30' region surrounding the 
bipolar system. All sources with excessive emission in the near- 
infrared are plotted in red. The exciting source of the bipolar 
system is presented as an asterisk, which shows prominent ex- 
cessive emission in the infrared and is in line with the protostel- 
lar nature of the system. 

Fig.E(fo) is the [3.6]-[4.5] vs. [5.8]-[8] map plotted with cat- 
aloged point sources from GLIMPSE I within the same region. 
The sample sources are shown as black dots. Sources with colors 
of 0.2 < [3.6] - [4.5] < 0.8 and 0.4 < [5.8] - [8] < 1 .5 are taken 
as Class II-type, while those wit h colors of [3 .61 - [4.5] > 0.4 and 
[5.8] -[8] ^ 0.4 as Class-I type dMegeath et aLOOO^lAlien et al.l 
12004 120071) . Presented as an asterisk, the exciting source of the 
bipolar system is confirmed as a Class I source. 

A sample of evolved and young stellar objects are plotted on 
a log[F60/F25] vs. log[F25/F12] diagram (Fig.|4]c). The central 
source of the bipolar system, presented as a blue asterisk, is lo- 
cated in a region dominated by young stellar objects, Hn regions 
and/or Herbig Ae/Be stars. The position of the exciting source 



1 http ://caravan. astro, wise .edu/protostars 



of the bipolar nebula provides additional evidence for the proto- 
stellar nature of the system. 

6. Conclusions 

We reported the discovery of an hourglass-shaped nebula at 
(RA.=18 A 14 m 24'.56, Dec.=-18°24'42".69 J2000) based on 
GLIMPSE data. The bipolar system is oriented in the north- 
east and south-west direction. A distance of 2.41 kpc is esti- 
mated based on galactic extinction models. A dust temperature 
of Td = 29.8 K and a luminosity of L wta \ = 5507 L are derived 
from the IRAS fluxes, which suggests the ongoing formation of 
a massive star or a quite massive stellar system in this region. 

Archived Bolocam 1.1 mm data, tracing cold dust, indicate 
the existence of an extensive structure at the waist of the bipolar 
system. The mass of the BGPS source is estimated to be around 
69 Mo, which is in line with the mass estimate of the envelope 
from the SED fitting (M env a 22 M Q ). This substantiates the 
detection of a massive envelope associated with the outflow sys- 
tem. 

The SED fitting results in a spectral index of 1.2 in the 
near- to far-infrared, which, along with the high infrared excess, 
strongly suggests the detection of a protostellar system. A stellar 
mass of 9.2 + 0.8 M Q and a luminosity of (4 + 1) x 10 3 L Q is 
obtained, which is consistent with the 5507 L Q estimated from 
the IRAS fluxes. 

This study substantiates the discovery of a massive protostel- 
lar system driving a bipolar outflow. Future investigations with 
higher spatial resolution may provide helpful information on the 
the kinematics, dynamics and chemistry of the inner and denser 
region, and contribute to our understanding of massive star for- 
mation in its early stages. 
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